Stochastic field theory for a Dirac particle propagating in gauge field disorder 
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Recent theoretical and numerical developments show analogies between quantum chromodynam- 
ics (QCD) and disordered systems in condensed matter physics. We study the spectral fluctuations 
of a Dirac particle propagating in a finite four dimensional box in the presence of gauge fields. 
We construct a model which combines Efetov's approach to disordered systems with the principles 
of chiral symmetry and QCD. To this end, the gauge fields are replaced with a stochastic white 
noise potential, the gauge field disorder. Effective supersymmetric non-linear cr-models are ob- 
tained. Spontaneous breaking of supersymmetry is found. We rigorously derive the equivalent of 
the Thouless energy in QCD. Connections to other low-energy effective theories, in particular the 
Nambu-Jona-Lasinio model and chiral perturbation theory, are found. 
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The propagation of an electron in a finite sample, such 
as a piece of wire, becomes diffusive due to multiple scat- 
tering at the impurities. The ensemble of the impurities 
is referred to as disorder. The diffusion constant of this 
process determines, together with the size of the sample, 
a universal energy scale, the Thouless energy. This, in 
turn, measured in units of the single particle mean level 
spacing yields the dimensionless conductance of the wire, 
see the review in Ref. [Q. Recently, it has been argued 
that the Gell-Mann-Oakes-Renner (GOR) relation has 
for QCD an analogous meaning &-0] , implying the exis- 
tence of a Thouless energy in QCD. Since the Thouless 
energy sets the scale within which the fluctuation prop- 
erties are fully of Wigner-Dyson type, it can be found by 
analyzing the spectral statistics. Indeed, this scale was 
identified in data of lattice gauge calculations gH . 

Our goal is the description of the diffusion process 
for QCD by a stochastic field theory. To this end, we 
merge Efetov's supersymmetric approach to disordered 
systems H with the principles of QCD and chiral sym- 
metry. This extends and complements the semi-classical 
reasoning of Refs. [0-0. We will not use formal analo- 
gies to the GOR relation as in Refs. [0-0, rather we will 
derive the Thouless energy as a natural result from our 
stochastic model. Anticipating our analytical findings, 
we draw an intuitive picture in Fig. |l|. Due to spon- 
taneous breaking of chiral symmetry, the gluonic gauge 
fields in the Yang-Mills action generate closed gluon and 
fermion loops in the QCD vacuum. We may view these 
vacuum fluctuations as "gauge field disorder" for the 
propagation of a constituent quark or, more generally, of 
a Dirac particle in a four dimensional box of finite volume 
V4. The multiple scattering at the vacuum fluctuations 
renders the motion of the quark diffusive. In this way, 
the pion decay constant and the constituent quark mass 
are generated. 

For our stochastic model, we introduce the Euclidean 
Dirac operator in four dimensions 
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where Ci, i — 1,2,3 are the Pauli matrices, CT4 — 
I2, and g is the gauge coupling. The complex fields 
u^i{x) describe the gauge field disorder. We choose 
them as stochastic white noise potentials with moments 
{Uf_,{x))u = and {u^{x)uy{y))u = ^5fj,^5^'^\x - y). The 
strength constant 7 will be determined later. Thus, the 
average of a quantity TZ[u] is given by 



{Tl[u])u = / d[u\R[u] cx]) ( / d'^xu^{x)u'^{x) 



(2) 



This replaces in our model the average over the non- 
Abelian gauge fields in the Yang-Mills action. 




FIG. 1. Propagation of a quark (thick solid line) through 
the QCD vacuum in a box of volume V4. The vacuum fluc- 
tuations involve quarks (thin solid lines) and gluons (twisted 
lines). The propagating quark is scattered at these vacuum 
fluctuations. 



Thus, we abandon gauge invariance, as can be seen 
from Eq. (||). We are free to do so because we are exclu- 
sively interested in stochastic features of the spectrum. 
The structure of the original gauge group enters in the 
model only through the assumption that the Dirac op- 
erator is proportional to the unit matrix TIn^ in color 
space. This is fully consistent with standard models of 
QCD !§. 

The eigenvalue equation iD[u\ipi = Ai[w]'0i defines for 
each realization of the disorder potential u{x) a spectrum 
{Ai[M]}- We wish to discuss the spectral correlation func- 
tions of k eigenvalues \p,p — 1 , . . . , fc 
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The energies are given proper imaginary increments to re- 
move the singularities on the real axis. The functions (H) 
can be obtained as derivatives of a generating function 
with respect to external sources Ji, . . . , Jfc, 
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where we have introduced A — diag(Ai, . . . , A^) ® I2 and 
J = diag(Ji, . . . , Jfc) (g) k with k — diag(— 1, -1-1). In the 
supersymmetry method, the generating or partition func- 
tion is expressed as a functional integral over the 4 • 2k- 
component and space time dependent superfields xi^) 
through 

Zfc(A + J) = (detgi?[u;A + J])„ 

d[x] exp (i f d^xx\x)D[u; A + J]x{x)) ) (5) 



with detg denoting the superdeterminant. Here, we have 
defined 



D[u; A + J] = (A + J) (g) I4 - l2fc <E) iD[u] 
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We suppress a Kronecker product with Ijv^ ■ The entries 
of A + J can be viewed as virtual valence quark masses. 
Due to the chiral structure of the Dirac operator the su- 
perfields xi^) can be decomposed in left and right handed 
parts Xl.,r{x) — Pl,rx{x), where Pl.r are the left and 
right handed projectors acting on the Dirac spinors. The 
disorder average according to Eq. (^ can now be per- 
formed and yields a field theory quartic in x{^)- Similar 
to the steps taken in Ref. ||l^, it can be decoupled via 
a Hubbard-Stratonovich transformation by introducing 
2k X 2fc complex supermatrix fields u{x) such that the 
fields x(a;) can be integrated out, 



Zk{\ + J) = 2^'=' I d[a]d[crt] exp {-F[a, ct^ A + J]) 



Both, (t(x) and its complex conjugate a'^{x) are indepen- 
dent variables. This reflects that xl{x) and X-r(^) are 
independent degrees of freedom in the chiral symmetric 
phase. The free energy reads 

F[CT,at;A + J] =tTg I d'^x(_^^^a{x)a^x) 

+trlogi:'[a,CTt;A + J]) . (8) 



The trace tr = trDiractrcoior in front of the logarithm is 
taken over both Dirac and color indices, the supertrace 
trg is defined in the space of supermatrices. The super- 
matrix analogue of the Dirac operator is given by 

D[a,a^;X + J] = a{x) ® Pr + a^x) ® Pl 

+ (A + J) ® I4 - l2fc ®{i^) , (9) 

and describes a Dirac particle coupled to the superma- 
trix fields 17(0;). An important observation from Eq. (H) 
is that (t(x) is coupled only to the right handed modes 
of the Dirac spinors and ct^(x) only to the left handed 
modes. This implies chiral symmetry, since <j{x) and 
a^ {x) are independent variables. 

In contrast to supersymmetry in standard high energy 
physics, neither x{^) nor a{x) directly represent physi- 
cal particles, such as quarks or mesons. Rather, these 
fields are the degrees of freedom in the stochastic model. 
Nevertheless, the functionals defined in Eqs. (^ and (H) 
posses certain symmetries which can be interpreted as 
the supersymmetric version of particle symmetries known 
from QCD. The Lagrangian in Eq. (^ is invariant under 
the chiral transformations 

Xl{x) — » ulXl{x) , ulul^ = ul^ul = l2fe 

XbXx) > UrXr{x) , UrUr^ = Ur^Ur = l2fe , (10) 

with ul.r. G SU(fc/fc). We notice that U(fc/fc) is likely 
to lead to axial anomalies which we do not discuss here. 
The corresponding supersymmetric version of chiral sym- 
metry is S\J L{k/k) g) S\J R{k/k). The generating fimc- 
tional (M) exhibits again chiral invariance 



(t{x) — > ULff{x)uR and a\x) 



t(x)4 , (11) 
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if the explicit symmetry breaking is turned off, i.e. A+J = 
0. This is not true for non-zero values of the virtual 
quark masses, A -I- J 7^ which explicitly breaks chiral 
symmetry. However, there is a remnant symmetry. Con- 
sider the case of k degenerate imaginary valence quark 
masses, Ai = . . . = Afc = A, i.e. A -I- J = Al2fc, where the 
sources are set to zero. Then the functional is invariant 
under the transformation 

cr{x) -^ ua{x)u* and a* [x) -^ ua* {x)u* , (12) 

with we S\]{k/k). 



We observe a close relation between the one-point 
function in the microscopic regime (near zero virtuality) 
and the two-point function in the bulk (far away from 
zero virtuality). In both cases, the symmetry breaking 
is controlled by a single parameter, although the sym- 
metry groups are different. In the bulk, the symmetry 
is broken by non-degenerate eigenvalues Ai and A2, and 
the symmetry breaking parameter is w = A2 — Ai . In the 
microscopic region, the symmetry is already broken by a 
non-zero eigenvalue A = Ai. 

In disordered systems, to make analytical progress, one 
considers the field theory via a saddle point approxima- 
tion in the limit of weak disorder, i.e. for a relatively low 
density of impurities. This corresponds to a Born approx- 
imation of the self-energy where diagrams with crossed 
interaction lines are neglected 0. Here, we proceed sim- 
ilarly by a saddle point approximation of the partition 
function (||) in the limit of weak gauge field disorder. 
This coincides with a 1/Nc expansion, which is a stan- 
dard asymptotic limit in many QCD models. These two 
limits coincide because for Nc ^ 00 and Ncg^ = const 
[ jll| , [l2[ the relevant diagrams contributing in leading or- 
der to the self-energy within our model are analogous to 
the ones in disordered systems. This shows the equiva- 
lence of the 1/Nc expansion and the weak disorder limit. 
From the saddle point condition SF[(7, u'l'; A + J] = 0, we 
find the self-consistent matrix equation in momentum 
representation 
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with H = (A + J) (Ki 14 and Eq = ctq «) Pr + crl® Pl, 
where cto smd crj are the fields at the saddle point. A 
cutoff Acut has to be introduced to regularize the diver- 
gent momentum integration. For scalar fields, the saddle 
point equation (Oh coincides with the gap equation of the 
Nambu-Jona-Lasinio (NJL) model p3[. The NJL model 
is a theory for chiral fermions coupled via a four-point 
interaction. The gap equation describes the self-energy 
of a quark calculated in the self-consistent Hartree ap- 
proximation obtained from a 1/A^c expansion H. For 
vanishing virtual quark masses, A -f J = 0, there ex- 
ists only a non-trivial solution, indicating chiral sym- 
metry breaking, if the coupling exceeds a critical value 
73^ > 2n'^ /NcAl^^.. For non-vanishing masses, however, 
one complex solution always exists, which we denote by 
M^. Since the spectral density per four-volume is given 
by v{\) = Im_Ri(A), see Eq. (0), the value of the con- 
stant 7 in Eq. (0) is now fixed. From Eq. (|l3) we find 
7 = ImM*/(7rg2^(A)).^ 

In contrast to the situation in disordered systems, we 
have to distinguish the microscopic and the bulk region. 
In the latter, the symmetries of the saddle point solution 



are still the same as given in Eq. (|lj). A solution of the 
saddle point equation is 



<7a = era 



ilmM^L 



(14) 



where L = 14(g) fc. On the other hand, in the microscopic 
region the chiral symmetry defined in Eqs.([lO|) to (11) is 
not preserved at the saddle point. The solutions read 

ao{x) ^ iM*Uo{x) and aoHx) ^ iM*ul{x) (15) 

where Uo{x) is 2 x 2 unitary supermatrix field and M* is 
the non-trivial scalar solution of the gap equation (|l^) 
for A + J = 0. We conclude that chiral symmetry is 
spontaneously broken such that 



SVR{k/k) (g> SVL{k/k) — > SV{k/k) 



(16) 



This is the supersymmetric analogue of spontaneous 
breaking of chiral symmetry in our stochastic model. 

By integrating out the quadratic fluctuations around 
the saddle points, we obtain, first, a theory for the two- 
point correlations in the bulk and, second, a theory for 
the one-point function in the microscopic region. In 
both cases, the zeroth order term is expanded up to lin- 
ear order in the symmetry breaking parameter, i.e. u in 
the bulk and A in the microscopic region, respectively. 
Moreover, we restrict ourselves to fields varying slowly 
in space-time, allowing for a gradient expansion of the 
quadratic fluctuations around the saddle points. In the 
bulk we obtain for the generating functional of the two- 
point function 

Z2{X,Lu,3) = J d^i{Q)c,cpi-F[Q;X,Lo,3]) (17) 

with 

F[Q: A, uj, J] = |i^(A)trg / d^x {V{\)d^Q{x)d^,Q{x) 

-2i{{u + ie)L + 3)Q{x)) , (18) 

where A = (Ai + A2)/2. The fields Q{x) are the Gold- 
stone modes of the saddle point manifold. They are in the 
coset manifold C/(l, l/2)/(t/(l/l) ® U{l/l)). Here, L is 
the symmetry breaking matrix and ie a proper imaginary 
increment to ensure convergence. This formally extends 
Efetov's result Q] to four space time dimensions. In par- 
ticular, the coset spaces coincide. Chiral symmetry is not 
present. This is in agreement with its explicit breaking 
by the virtual current mass i\. The diffusion coefficient 
I?(A), however, a main object of our interest, must be dif- 
ferent from disordered systems because we start out from 
the Dirac equation. It is given as the second moment 

tr / d'^y{x - yfS{x, y;X + ilmM^)S{y, x;X- ilmM^) 



2TTiy{X)V{X) 
(ImAf;)2 



(19) 



where S{x,y;X — iluiM^) is the Green function to the 
mean field Dirac operator A + iluiM^ — i ^. 

From the competition between the kinetic and the sym- 
metry breaking term one can infer the Thouless energy 
in the bulk. In a finite volume V4, the longest wavelength 
which fits inside the Euclidean box is given by the inverse 
of the linear extension of the box, i.e. 1/(V4)^/^. Thus, 
the kinetic energy cannot be smaller than I?(A)/ y^- For 
smaller energies, the wavelength exceeds the extension of 
the box and does not resolve details inside it. In that 
case, the fields become constant and the functional is 
dominated by the constant modes. Hence, only the sym- 
metry breaking term gives a contribution. Then the level 
statistics is of Wigner-Dyson type. If the energy exceeds 
this scale, the kinetic term becomes more and more im- 
portant, leading to deviations from the Wigner-Dyson 
statistics. Thus, the equivalent of the Thouless energy in 
the bulk is given by 

V bulk 
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V{X)/y/Vi 
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In the chiral limit, A — > 0, the diffusion coefficient (|I 
in momentum representation reduces to 

8NciM*y f d'^p 1 



V{Q) 



IW)I 



(2^)4 (p2 + (Af *)2)2 



(21) 



where (V'V') is the chiral condensate, obtained from the 
Banks-Casher relation, (■!/'■!/') = — Trlimy^^oo ^(0). For 
this derivation, we used only the basic assumption (g) 
and chiral symmetry. Low energy theorems of QCD im- 
ply that the integral in Eq. ( pl| ) turns out to be propor- 
tional to the square of the pion decay constant /^, 



fl^AN,{M*f 



d'^p 



1 



(22) 



for a detailed discussion see Ref. H 
information, we arrive at 

P(0)=2/2/|(ViV)| 



(2^)4 (p2 + (M*)2)2 ' 

With this additional 



(23) 

In Refs. P-B, this formula was obtained from the GOR 
relation by using formal analogies and semi-classical con- 
siderations. Here, we gave a rigorous derivation of the dif- 



fusion constant (21) from our stochastic model. Finally, 



we obtain for the generating functional of the one-point 
function in the microscopic regime 



(24) 



Zi{X, J)= dfi{U) exp i~F[U; X + J] 
with 
i^[C/;A + J]-|(V:;^)|trg /"d*x( 



■iw>r^ 



d^uix)df,uHx) 



-'-{X + Jk){U{x) + U\x))) . (25) 

The supermatrices U{x) are 2x2 unitary. Again, the 
functional is dominated by the constant modes for ener- 
gies smaller than f^/{\{^^)\y/V4). We conclude that the 



equivalent of a Thouless energy in the microscopic region 
is given by 



A" 



ViO)/^V^^2f^/{\{H^)WVi) 



(26) 



We note that the microscopic functional reduces, up to 
irrelevant numerical factors, to the partially quenched 
chiral perturbation theory (pqCPT). This is consistent 
with the considerations of Refs. |14| , p^ , were pqCPT was 
used to identify the equivalent of the Thouless energy in 
the Dirac spectrum. 

In summary, we presented a stochastic field theory by 
merging Efetov's approach to disordered systems with 
the standard principles of QCD and chiral symmetry. 
Due to the latter, we obtain different theories in the mi- 
croscopic and the bulk region which connect to the NJL 
model and pqCPT. The Thouless energy is a natural con- 
sequence of our model. In conclusion, we wish to make a 
proposal: since we now have a detailed and, importantly, 
quantitative description of the spectral fluctuations in 
QCD, it would be worthwhile to use this insight as an 
input for lattice QCD, either by constructing effective 
models or by modifying the algorithms accordingly. It 
would be highly desirable if this made lattice QCD more 
efficient. 
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